Convergence of the chiral expansion for the nucleon mass 
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A number of papers recently have used fourth-order chiral perturbation theory to extrapolate 
lattice data for the nucleon mass; the process seems surprisingly successful even for large pion 
masses. This paper shows that the inclusion of the fifth-order terms spoils the agreement. 
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Over the last few years there has been an explosion of activity in the field of chiral extrapolations of lattice QCD data. 
However for many quantities of interest unqucnched calculations have typically only been performed at relatively high 
quark masses, so that the pion mass is 500 MeV or more. There are real questions about the convergence of chiral 
expansions in this region. 

One quantity for which the chiral expansion has seemed to work surprisingly well is the nucleon mass. Various groups 
have looked at the 0(p 4 ) (technically N 3 LO) expansion of the nucleon mass in heavy baryon chiral perturbation 
theory (HBxPT) and found good agreement with data up to pion masses of 800 MeV or more. Of course the results 
are quite sensitive to the input parameters which include some rather poorly known low-energy constants, but when 
these are used as fit parameters the results agree well with other determinations. 

Some years ago the 0(p 5 ) contribution in HB^PT was calculated 0. At this order two-loop graphs enter, but it 
turns out that there are no irreducible two-loop contributions and all the reducible graphs can be absorbed into the 
physical pion mass and the pion- nucleon coupling constant in the third-order contribution. Only a relativistic (1/M N ) 
correction is left, and it is tiny: just 0.3% of the third-order term for = 139 MeV. 

The details of this calculation have not so far been used by any group exploring the chiral expansion of the nucleon 
mass. Refs. 0, include the relativistic correction (which can also be obtained from a fully relativistic one- 
loop calculation), while Ref. Q uses its smallness to support the view that 0(p 5 ) contributions will not spoil the 
convergence. However this is not correct. So long as one is interested in comparing the predictions of xPT to 
experiment in the real world, it is both legitimate and necessary recast the results of calculations, obtained initially in 
terms of bare or chiral limit parameters, in terms of their physical values. But if one wishes to explore the variation 
with quark mass in lattice calculations, the expressions must be left in terms of bare parameters. (If these are not 
well-known they may in practice be given the same numerical value as their real-world equivalents, but that is a 
different matter.) 

Ref. p] was not written with lattice extrapolations in mind, and so the appropriate expression was not given explic- 
itly. Furthermore the calculation was performed using the Lagrangian and notation of Ref. @. This non-minimal 
Lagrangian contains more parameters (31) than there are independent LECs (23), so the extra ones are taken to 
vanish at a particular renormalisation scale - the physical pion mass. This is an unfortunate choice if the aim is to 
make explicit the running with (lattice) pion mass. The aim of the current paper is to present the 0(p ) result in 
an appropriate form for comparison with lattice results, and to explore its consequences for the convergence of the 
chiral expansion. We should remark that Beane has also opined that fifth-order corrections could be significant, using 
a rather different argument 0. 

The nucleon mass in the chiral limit is denoted by M. The next four terms in the chiral expansion for non- vanishing 
quark mass, and hence bare pion mass m, are as follows: 
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The parameters g, f and m are the leading terms in the chiral expansions of g A , f v and m w ; m 2 is proportional to the 
quark mass. The physical values of these parameters are g A hys = 1.267, /P hys = 92.4 MeV and mP hys = 139.6 MeV; 
in what follows the bare parameters g and / will be set to the physical values as the chiral-limit values are not well 
known. The second-order LECs Cj are more-or-less well determined from pion-nucleon and nucleon-nucleon scattering. 
Since this paper follows most closely the methods of Refs. [3j,|4| we use the same values as there: c\ is generally fitted 
but has been determined to be in the range —1.1 to —0.4 GeV -1 ; C2 = 3.2 GeV -1 and C3 = —3.4 GeV -1 ; the exact 
values of the latter two do not matter because they only occur along-side higher-order LECs which will be fitted. 

The finite third-order LEC dig is fitted to the Goldberger-Trieman discrepancy, and using g^ NN = 13.18 gives d\% = 
-0.78 ± 0.27 GeV -2 [fj. With this as input, the value of d r 16 (mP hys ) is extracted from the reaction ttN -> ttttN in 
Ref. 0; however the same authors subsequently use a substantially revised value, and Ref. quotes their 1 a range 
to be -2.61 to -0.17 GeV~ 2 . 

The notation for the mcsonic LECs is if (A) = (7;/32tt 2 ) (I, + 2 lnm^/A), with 7 3 = 2.9 ±2.5, 1 4 = 4.4±0.3, 73 = -5 
and 74 = 2@. 

By e(A) we denote a combination of LECs: e(A) = 16eg 8 (A) + 2e^ 15 (A) + 2e^ 16 (A) — 32cid2 S (\), where d{ and are 
third- and fourth-order LECs respectively. 

The fourth-order calculation was first given in Ref. [l0( (where fourth-order LECs were denoted bi). The fifth- 
order calculation is from Ref. with d^sM restored; it was omitted there because it vanishes at A = ?7iP hys . An 
unconventional notation for the mesonic LECs was also used in that paper. 

If lattice results presented the variation of M N with the quark mass, the expression above would be the relevant one. 
Since, however, the usual presentations plot M N against m n , both being lattice measurements, the correct expression 
uses the NLO — 0(p 4 ) — running of the pion mass to give 
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with e'(A) = e(A) + 8cil^(\) / f 2 . At fourth order it make very little difference which expression one uses. 

The LECs that enter can, in principle, all be extracted from irw, 7rN and NN scattering data and the physical nucleon 
mass. As indicated above, most already have been, with greater or lesser uncertainties. However the fourth-order 
LEC e'(A) is very poorly constrained and is generally fitted to lattice data. Bernard et al recognised that with one 
free parameter the fourth-order expression could fit lattice data out to m v of about 600 MeV [ll|; Procura et al 
showed that the best fit to the lattice data up to 600 MeV was actually a superb fit right up to much higher masses 
(though without much predictive power when one varied the fit parameters within acceptable limits of y 2 ) 4]. From 
such results it has been concluded that the chiral expansion may be valid much further than naive estimates might 
suggest. However until now the full fifth-order terms have not been included. 

In Fig. ^ we show the analogues of the fit of Procura et al Q| at third, fourth and fifth order ( Eq.El ). The lattice data 
is from the following collaborations: UKQCD and QCDSF 0, CP-PACS HG3, JLQCD [|[l|Und Orth, Lippert 
and Schilling 0] ; the selection of points at large enough lattice volume is made by Musch in Ref. |l5j ; see also Ref. Q . 

Like Procura et al, at fourth order we fit M and e'(m^ hys ) to the lattice data below 600 MeV, choosing c% so that the 
physical M N is reproduced. At third order in HBxPT only M is free. At fifth order we choose to fit c^ 6 (mP hys ) since 
it is poorly determined otherwise. However treating M (or equivalently c\) as a fit parameter gives the unacceptable 
solution M > M N , so instead we fix c\ at the value emerging from the fourth-order fit. 

We see that at third order an acceptable fit cannot be obtained: % 2 = 31.6 for 3 d.o.f. However the value C\ — 
— 1.14 GeV -1 (M = 865 MeV) is quite reasonable. At fourth order the fit is extremely good, with x 2 = 0.20 for 2 
d.o.f.. Also, the value c\ = —0.84 GeV -1 is within the range obtained by other means and e'(mf ys ) = —0.52 GeV -3 
is natural in size. 1 At fifth order the fit is equally acceptable: % 2 = 0.54 for 2 d.o.f. However the fitted values are well 
out of line with what one would expect: d^ 6 (mP hys ) = 4.654 GeV~ 2 and e'(mf ys ) = —30.5 GeV -3 . Furthermore in 
sharp contrast to the fourth-order result, the fifth-order curve, though good for the points it was fitted to, breaks down 
at higher values of m^. Allowing the fit to vary within an acceptable range of y 2 does not change this behaviour. 



1 This is equivalent to e'(l GeV) = — 3.67GeV 3 . The slight deviation from the values found by Procura et al is due to the fact that they 
did not include the corrections for the running of and also to their inclusion of the "naive" fifth-order term. 
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FIG. 1: Best fits to the lattice data (constraint to pass through the physical point) for M N (in GeV) versus m 2 (in GeV 2 ) 
below = 600 MeV (unshaded region) at third- (blue, long dashes), fourth- (red, short dashes) and fifth order (black, solid). 
See text for the sources of the data. (Note that the apparent agreement of the fourth- and fifth-order curves at low ra 2 masks 
very different fitted values of the LEC e'.) 

Whereas the fourth-order fit results in distinctly small coefficients, M N = 0.887 + 3.35mJ — 5.61mJ + 3.96m 4 . (in 
GeV), the fifth-order fit does not: M N = 0.885 + 3.20m; 2 . - 5.61tt4 + 33.9m 4 - 108. 5m^. (The logs are evaluated at 
to^ = 0.2 MeV but do not change the qualitative picture). In the latter the large size of the fifth-order term is not 
an artifact of the unrealistic value of d\ 6 (m^ hys ); with this set to zero the coefficient is still 56. (Interestingly though, 
the order of magnitude is not so very different from that found (« —40 GeV -4 ) in most of the finite-range regularised 
expressions used to extrapolate over a very large range of pion masses by Leinweber et al ^(| . 

As already mentioned, the bulk of the fifth-order term comes from the running with m^ of the pion-nucleon coupling 
constant, and in particular of g A . But it is well known that NLO xPT does not do a good job of reproducing lattice 
data on g A ; the best fit to date has included an explicit A |17|. which is absent here. As an experiment, we added 
an m% term to the expression for the nucleon mass as if it came from the m^ term in g A , with a coefficient chosen 
to give a reasonable fit to the lattice data presented in Ref. 0. 2 Alternatively we simply dropped all terms in the 
running of M N originating from the running of g A (since the latter is apparently quite weak). The results are shown in 
Fig. |21 Although the resulting curves are quite different from the standard one, they cannot be said to be much more 
successful at extrapolating beyond the fit region, nor are the corresponding extracted LECs notably more acceptable. 



2 This is achieved with (ijgfm* 5 ) = 1.24 GeV -2 and (in g A only) c 3 - 2c 4 = 0.62 GeV -1 . The third-order expansion of g A was first 
given in Ref. [H. 
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FIG. 2: Best fits to the data below = 600 MeV at fifth order (black, solid), without allowing g A to run (magenta, long 
dashes), and with a 6th-order term designed to fit the lattice running of g A (green, short dashes). See text for more details. 

Ultimately, we do not consider the results of this paper to be a surprise. What was surprising was the original finding 
that HBCPT appeared to work not only out to 600 MeV but even beyond. Our results imply that lattice data will 
be required at substantially lower pion masses before chiral extrapolations will be reliable. 
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